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PREFACE 

This little book was written at the instance of Messrs. Adam 
Hilger, and, in accordance with their desire, it contains just what 
is required for the purpose of reading and handling my Simplified 
Method of Tracing Rays, etc. (Longmans, Green & Co., London, 
1918). With this practical aim in view, all critical subtleties have 
been purposely avoided. In fact, it is scarcely more than a synop- 
tical presentation of the elements of Vector Algebra covering the 
needs of those engaged in geometrical optics. At the same time, 
however, it is hoped that this booklet will serve a more general 
purpose, viz. to provide everybody unacquainted with the subject 
with an easy introduction to the use of Vector Algebra. 

It is scarcely necessary to explain that the deductions given 
in this book are based on Euclid's axioms, notably with the 
inclusion of his postulate of parallels — upon which the equality of 
vectors is most essentially based. Those readers who are desirous 
of seeing how the formal rules here given can be generalized so as 
to be valid independently of the axioms of congruence and of 
parallels, may consult the author's Projective Vector Algebra (Bell 
& Sons, 1919), and a sequel to it published in Phil. Mag. for July, 
1919, pp. 115-143. It is, however, advisable for the student to 
become first thoroughly familiar with the euclidean vector algebra 
as here presented. 

I take the opportunity of expressing my sincere thanks to Messrs. 
Hilger for enabling me to make this further contribution towards 
the promotion of the more general use of this powerful and 
convenient language of vectors, and to the Publishers for the care 
they have bestowed upon this little book. L. S. 

London, August, 1919. 
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ELEMENTS OF VECTOR ALGEBRA 

1. Vectors defined. Whereas common algebraic magnitudes, 
such as the number of inhabitants of a village, or the mass of a 
body, or the energy stored in an accumulator, having nothing to 
do with direction, are called scalars, any magnitude such as a 
displacement, a velocity or an acceleration, which has size as well 
as direction in space, is called a vector. The visual, or tangible, 
representative of any vector whatever is a segment of a straight 
line of some length, representing the vector's size, and of some 
definite direction in space, together with its sense (say, from a 
point M towards a point N), giving the direction of the vector. 

Vectors will be printed in Clarendon, thus 

A, B, etc., or n, r, s, etc., 

and their sizes, regardless of direction, or their tensors (as they 
are called) will be denoted by the same letters in Italics. Thus, 
A will be the tensor of A ; B, n will be the tensors of B, n, and 
so on. 

Returning once more to the above definition, we may as well 
say that any vector A = OE is given by the ordered couple or pair 
of points, the origin and E the end-point of the vector ; the 
tensor, called also the absolute value, of the vector being 
the mutual distance of and E. In short symbols, and using 
the familiar bar for the distance, 

A = O^E, A = OE. 

The tensor of a vector is thus an ordinary, absolute or essentially 
positive number. 

A vector whose tensor is (in a conventionally fixed scale) equal 
to unity, is termed an unit vector. Thus, if r = I, the corresponding 
r will be a unit vector. It will be understood that the denomi- 
nation of A is that of A. That is to say, if A is, for instance, the 

S.V.A. a 
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displacement of a particle, A will mean so many centimetres ; 
and if A represents a velocity, A will be a number of cm. per second, 
and so on. 

As far as will be possible we shall reserve small (in distinction 
from capital) Clarendon letters for unit vectors. Thus, if the 
contrary is not expressly stated, a, b, etc., will stand for unit 
vectors, so that a — i, b = i, etc. 

In MS. work the reader will, at least in the beginning of his 
vector career, find it useful to underline all his vectors once or 
twice. Or he may write them thicker, and imitating somehow the 
printer's type. Then, everyone will soon find out his most agreeable 
manner of writing. 

2. Equality of vectors defined. We have just seen that the 
two essential features of a vector are : its size or tensor, and its 
direction in space. "~~- 

In some branches of physico-mathematics it is important to 
consider the position of the vectors in question (besides their sizes 
and directions), i.e. to localize their origins, either by fixing the 
origin of each vector altogether or by allowing it only to move 
freely in its own line. Such vectors are usually called " localized " 
vectors. In a vast class of investigations, however, the position 
of these directed magnitudes is of no avail, and it is then obviously 
convenient not to include . position among the determining char- 
acteristics of a vector. Such vectors, in distinction from localized 
ones, are called free vectors. These and these only will here 
occupy our attention. The adjective will be dropped, however, 
and the beings in question will be called shortly vectors. With 
this understanding, the definition of their equality may be put 
thus : 

By saying that two vectors, A and B, are equal to one another, 
and by writing 

A=B or B=A, 

we mean that their tensors are equal, A=B, and that they have 
the same direction or, in other words, that the straight segments 
representing these two vectors have the same length and are con- 
currently parallel to one another. In short symbols, A = B means 
as much as 

A=B and AttB. 
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Thus, if a pair of points, 0, E, represents a vector A = -* E, 
the oo 3 point pairs 0', E' or straight segments O'E' of equal length 
with and concurrently parallel to OE are all equal to A, no matter 
where their origins are situated. Notice that through every point 
0' of euclidean space there is one and only one parallel to OE, 
so that from every space point 0' as origin one and only one vector 
can be drawn which is equal to the given A. Of course, the 
laying off, from 0', of the length 0'E' = OE implies the use of 
some " rigid transferer," such as a pair of compasses. 

Equivalently, we may say that the rigid translation (parallel 
shifting) of a given vector is irrelevant, or does not change the 
vector. Provided it is not being rotated, stretched or contracted, 
we can, by the accepted definition, " transfer " it to any place we 
like best. 

Two vectors A, B drawn from the same origin are termed 
coinitial. By what has just been said, any two vectors can be 
made coinitial, by shifting one of them or both parallel to them- 
selves. If A=B, then making them coinitial, fuses them into 
one straight segment. If only A=B (equal tensors only), then 
making the vectors A, B coinitial will still leave a certain non- 
vanishing angle, or direction difference, between them, sufficient by 
itself to declare the two vectors as being different from one another. 

We will say that two or more vectors form a chain if the end- 
point of one serves as the origin for the other, and so on. As 
before, any two vectors A, B can be linked up into a chain, to 
wit in two manners : end-point of A coinciding with origin of B 
(or A preceding B), or vice versa. 

This licence will be seen to be of capital importance for the 
vector sum to be defined presently, inasmuch as it will confer 
upon that sum the extremely convenient property of commuta- 
tivity. It will, therefore, be important to keep these latter, 
apparently trivial remarks well in mind. 

3. Addition of vectors. Let A and B be any two vectors, 
drawn anywhere. Shift B so as to bring its origin to coincidence 
with the end-point of A, as shown in Fig. 1. The vectors being 
thus linked up into a chain we call sum of A and B and denote by 

S=A+B 
a third vector S which runs from the beginning to the end of the 
chain, i.e. from the origin of A to the end-point of B. 
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This is the definition of the vector sum. The operation, vector 
addition, thus defined has the so-called group-property, that is to 
say, being performed on vectors it gives again a vector, in much 
the same way as five apples added to three apples give again a 
certain number of apples. 




Fig. i. 

The above vectorial expression will be read : B added to A. 
But we might as well have linked the two given vectors so that 
the end-point of B=B' falls into the origin of A, as shown in the 
lower part of Fig. I. Then their sum, say S', would — according 
to the definition — be 

S' = B+A, 

which reads : A added to B. The natural question arises : What 
is this new vector S' ? Is it equal to S ? 

The answer is in the affirmative. For, by construction, B' is 
parallel to B and B' = B, so that OaB and aOy are congruent 
triangles, and S' = S. At the same time the angles jSand y are 
equal to one another, so that, a/3 being parallel to yO, so are also 
Of$ and ya, or S MS'. Therefore, by Section 2, S' = S, what was 
to be proved. 

Thus we have 

A+B=B+A, (i) 

the commutative property of vector addition. The order of the 
addends, in the vector chain, is irrelevant for their sum. 

Again, we might have shifted B to the position B" (Fig. i), 
retaining also the previous B = a->/3 and constructing A r — 5-*-/? = A. 
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Then, 0a/3& being a parallelogram and S = -> /3 one of its diagonals, 
we should have the following construction of the sum of two 
coinitial vectors A, B (Fig. 2) : 

Through a, the end-point of A, draw a parallel to B, and through 
/?, the end-point of B, draw a parallel to A. Then y, the cross of 




these parallels, will be the end-point of the required vector sum 
A+B or B+ A, and the common origin of the two addends will 
be the origin of their sum 

S = 0->y=A+B=B+A. (2) 

This is known as the parallelogram construction of a vector sum. 

We might have started from it as a sum definition. It has the 
advantage of being immediately symmetrical with respect to the 
two addends. At any rate we see that the chain and the parallelo- 
gram constructions are (in virtue of Euclid) wholly equivalent to 
one another. 

Thus far the case of two vector addends. Now, the sum of 
these being again a vector, S = A+B, we can add to S any third 
vector C, thus obtaining 

S + C = (A+B)+C=C + (A+B), 

the latter by the commutative property. Similarly for the sum 
of four and more vectors. Again, linking up the vector addends 
A, B, C into a chain, we see without difficulty (Fig. 3) that 

(A+B)+C=A + (B+C), (3) 

the result being in both cases the same vector, viz. that drawn 
from the beginning to the end of the chain. The same property 
holds for the sum of any number of vectors. The brackets become 
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superfluous, and either of the above expressions can simply be 
written 

A+B+C 
or B +A +C, and so on. 




The addition of vectors is thus seen to be associative as well as 
commutative, exactly as the ordinary algebraic addition of scalars. 

If by any appropriate parallel shifting of any number of given 
vectors, say A, B, C, D, they can be linked up, as in Fig. 4, into 




a closed chain (or a polygon), plane or not, then the sum of these 
vectors is a nil vector or simply nil, 

S=A+B+C+D = 0, 

and therefore also A + C+B+D=0, etc. It is scarcely necessary 
to say that a vector is nil or zero, S=0, if 

S=o, 

that is, if its tensor vanishes, and conversely ; or, in other words, 
if its end-point and origin coincide, such precisely being the case 
of our closed chain. 

The vector sum, which shares with the ordinary algebraic sum 
the two capital properties of commutativity and associativity, 
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contains the algebraic sum as a particular sub-case, to wit, when 
the vector addends are all parallel to one another. For, such 
being the case, they can always be brought into one line or made 
collinear. Parallel vectors, no matter what their tensors, are 
therefore called also collinear vectors. Now, if A, B are collinear 
vectors, the tensor of their sum is 

A±B, 

according as A, B are of equal or of opposite senses. 

The tensor of a sum of vectors, as S=A+B, can conveniently 
be denoted by 

S=|A+B|, 

as is usual for the absolute value of ordinary algebraic magnitudes. 
Thus we shall have, for collinear vectors, 

|A+B| = |^±B|. 

But, it will be well kept in mind that in general, for non-collinear 
addends, 

|A+B|,fe|^±£|, 

since |A+B| is the length of the third side of a triangle whose 
two other sides are A and B. 

By what has just been said, the sum of two equal vectors, which 
is written 

A+A or 2A, 

is a vector coinciding with A in direction and having 2A for its 
tensor. Similarly for 3A, 4A, and so on. 
Again, if B be such a vector that 

2B=A, 

we shall write 

B = JA, 

and similar meanings will be attached to JA, JA, etc. In this 
manner, and applying in the case of irrational factors the well- 
known limit-reasoning, we easily obtain the meaning of the 
expression 

wA, 

where n is any positive scalar number, integral, fractional or 
irrational. We can say shortly that wA is the vector A stretched 
in the ratio n : l. If n is negative, then (as justified in Section 4, 
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infra) mA will be the vector A stretched in the ratio \n\: I and 
then inverted in its sense, or first inverted and then stretched. 

In particular, if a is an unit vector, the " unit of A," as we 
have said before, we shall obviously have 

A = ^a. (4) 

Here A, the tensor of A,' is an ordinary positive number. 

Let a, b be any two non-collinear unit vectors. (Imagine them 
shifted so as to be coinitial.) Then any vector R contained in or 
parallel to the plane a, b can obviously be expressed by 

R=xa + yb, (5) 

where x, y are some scalar numbers. For the plain meaning of 
this assertion is that, starting from 0, the origin of a, b, any other 
point of the plane a, b can be reached by making a number [x) 




Fig. 5. 

of steps a and then a number (y) of steps b (or first yb and then 
xh). If both x and y are positive, then, with as origin, R will 
lie in the region I of the plane a, b (Fig. 5) ; if x < O, y > o, it 
will fall into II ; if x < o, y < O, into III, and finally, if x > 0, 
y < o, into the region IV. 

The scalars x, y in (5) are called the components of R along 
a, b as axes. 

Similarly, if a, b, c be any three non-coplanar vectors * which 
we may again take as unit vectors, then any vector whatever can 
be expressed in the form 

R=xa + yb+sc. (6) 

* I.e. such as cannot be made coplanar by parallel translations. 



ADDITION OF VECTORS 9 

The scalars x, y, z are called the components of R taken along 
a, b, C as axes. These axes may be chosen at our will (if we wish 
at all to split our vectors R into components), either perpendicularly 
or obliquely to one another, the only condition for covering all 
possible vectors (R) being that a, b, c should not be coplanar. 

The three vectors a, b, c or, as we will say, the reference system 
a, b, C, being fixed conventionally, we see from (6) that any vector 
is fully determined by three scalar data x, y, z and not less than 
three. The same thing is obvious from formula (4), according to 
which any vector R can be represented by 

Rt. 

In fact R is one scalar number, and r (a direction) implies two 
more scalar data, for instance two angles, which makes in all three 
independent scalar data as above. 

In so-called polar coordinates, for instance, we have (Fig. 6) 

R = i?{[icos <£+j sin <£] sin 0+kcos 6}, (7) 

where i, j, k are mutually perpendicular unit vectors, the pole- 
distance or co-latitude being counted from k and the longitude 
from i towards j, if i, j, k be a right-handed system. The particular 




Fio. 6. 



form (7) will often be found useful in passing from vector to scalar 
formulae, especially in optical computations. The unit of R, i.e. 
any unit vector r, will be, in polar coordinates, 

r = [icos <j> + j sin $1 sin + k cos 6. (7a) 



io ELEMENTS OF VECTOR ALGEBRA 

After what has been said, it is scarcely necessary to explain 
that every vector equation is equivalent to three scalar ones. For 
to give R is equivalent to giving, say, its three rectangular or 
" cartesian " components, or, as in (7), the polar coordinates 
R, 0, 4> of the end-point of R, with as origin. Thus A =B means 
as much as A 1 = B 1 and A 2 = B 2 and A 3 = B 3 , if the suffixes I, 2, 3 
are used for the components of the vectors along i, j, k, or as 
much as 

if the suffixes A>Ji are used to distinguish the polar coordinates of 
the end-point of A from those of the end-point of B. 

But it must henceforth be urged that any such splitting of a 
vector should be avoided as much as is possible in the course of a 
vector investigation of any kind. For the utility of vector method 
lies precisely therein that it enables us to treat vectors as wholes 
instead of the triad of " components " of each of them. 

4. Subtraction of vectors. This will require but a few remarks. 
In fact, as in common algebra, the difference of two vectors A 
and B, to be denoted by 

A-B, 

may be defined as such a vector C, which added to B gives A. 
In symbols, we say that 
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if B+C = 

From this definition we see at once that if A, B are made coinitial 
(Fig. 7), the vector A-B runs from the end-point of B to the 
end-point of A. From the same figure, and by what was explained 




previously, we see that A +B and A -B are represented by the two 
diagonals of the parallelogram constructed upon A, B. 
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Apply the above definition (8) to the particular case A=0; 
then 

C = o-B= -B, 
and B + C = O ; therefore, 

B + (-B)=o. (9) 

This settles the meaning of the vector denoted by -B ; it is the 
vector which runs from the end-point towards the origin of B, or 
the reverse of B. This also justifies the interpretation given 
before to a negative scalar factor of a vector. Henceforth, for 
any A, B, 

A-B 

wiil stand for the same vector as 

A + (-B). 

The above remarks complete the meaning of 

nA, 

where A is a vector and n any real scalar, positive, nil or negative. 
The concept of such a product of a vector by any scalar n does not 
contain, in fact, anything besides the previous concept of vector 
sum or difference. It is derived from their special case, viz. 
relating to collinear vectors. 

To say it once more, nA is simply the vector A stretched in the 
ratio \n\ : I and, if n < O, turned through l8o° (in any plane 
passing through A). 

Finally, as the reader himself will easily prove, for any A, B, 
and any scalar factor n, 

n (A+B)=wA+wB. (10) 

Similarly for three or more vector addends. This settles all 
questions concerning the multiplication (or division) of a vector 
expression by any scalar number. 

5. Scalar product of two vectors. We now come to a new 
concept, transcending that of vector addition which hitherto has 
occupied us. The " scalar product " of two vectors A, B which 
will be denoted by AB is, first of all, not a vector but a scalar. 
(Thus the scalar multiplication of vectors does not respect the 
group requirement ; it yields a result not contained in the class 
of operands : it takes two vectors and constructs out of them 
something which is utterly deprived of direction. None the less 
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it is a very useful operation.) The value of this scalar is, by 
definition, proportional to the tensors of both the factors and to 
the cosine of the angle (A, B) included between them. In short, 
the definition of the scalar product is 

AB = ^5cos (A, B). (i i) 

This can also be read : AB is the projection of A upon B multi- 
plied by B, or the projection of B upon A multiplied by A. 
Since AB = BA, for A, B are common numbers, and 

cos (A, B)=cos(B, A), 

we see at once from the very definition (il) that 

AB=BA, (12) 

the commutative property. 

According as the angle (A, B) is <^ or > - (but <— J, the 
product AB is positive or negative ; for A, B are themselves 
essentially positive. And if (A, B) = - or A J_ B, then 

AB=o, 

no matter what the (finite) tensors of A, B. In this case the 
operation (scalar multiplication) deprives the material operated 
upon not only of direction but of size. It annihilates it. 

Conversely, if of two vectors A and B we know only that AB =0, 
then the only conclusion we can draw from it is that 

AJ_B, 

but by no means that one of the factors vanishes, unless we happen 
to know beforehand that the two vectors cannot be perpendicular. 
It is of prime importance to keep this well in mind : 

AB =0 means in general only as much as A J_ B. 

The scalar product AB contains the ordinary algebraic product 
as a special case, to wit, when A, B are collinear vectors. For if 
such be the case, we have cos (A, B) = + I, and therefore, 

AB=±^5, (13) 

according as A and B have the same or opposite directions. 
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Since the tensor of the vector mh is mA, we see at once from 
(11) that 

raAwB = mnAB. 

Thus, for example, if a, b be the units of A, B, we have 

AB = ^Bab, 

where, again by the definition (11), 

ab = cos(a, b), (14) 

valid for any pair of unit vectors a, b. Thus, for instance, if a, b 

make with one another the angle of 45° we have ab=— 7=, and if 

V2 

(a, b)=90°, ab=0. For the three normal unit vectors i, j, k used 

above we have ij =jk=ki = 0. 

As a sub-case of (13) we have the scalar square of a vector, or 

better, its autoproduct, 

AA or A i = A i , 

and if a be a unit vector, 

a 2 = a 2 = 1. 
Thus, i 2 =j 2 =k 2 =i. 

Again, if R is any vector whatever and n a unit vector, Rn is 
the (scalar) component of R along n, or the orthogonal projection 
of R upon n as axis, 

Rn = i? cos (R, n). 

By what has been said we see that if A, B be rigidly linked 
together and thus moved about in space in any arbitrary manner 
whatever (spun round, etc.), the value of the product AB is not 
changed. It is thus an invariant of the pair of vectors with 
respect to their common rigid motion. In fact, AB depends only 
on the tensors of A, B and on their relative direction, i.e. the angle 
(A, B). 

By the fusion of A, B into AB all directional properties of the 
factors are gone. The result has nothing more to do with direction 
in space ; it is an ordinary scalar, like the tensor of each of the 
two vectorial factors. Thus, if C be a third vector, 

(AB)C or C(AB) 

will simply mean the vector C magnified (stretched) AB times, 
assuming, that is, that AB is a dimensionless or pure number ; if 
AB is an area and C, say, a displacement, then (AB) C, the tensor 
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of (AB)C, is a volume, of course, and so on. If D is a fourth 
vector, 

(AB)(CD) 

will again be a scalar, and so on. The brackets are here used as 
separators. They are, of course, indispensable in such and similar 
cases. For, to take only three factors, ABC would, in general, be 
ambiguous, since 

(AB)C is a vector along C, 
while 

A(BC) is a vector along A, 

and thus entirely different from the former. Instead of brackets 
dots may conveniently be used as separators, thus 

(AB)C=AB.C, 
(AB)(CD)=AB.CD, 

and so forth. The reader will soon find that this need of precaution 
gives rise to no serious inconvenience. 

The scalar product AB is commutative owing to the symmetry 
of its very definition with respect to A, B. In this it resembles 
the ordinary product. But, what is most important, it has also 
the distributive property, viz. for any A, B, C, 

A(B+C)=AB+AC. (15) 

For, by the definition, A(B+C) or (B+C)A is the projection 
of the vector B+C upon A multiplied by A. But the projection 
of the sum of two (or more) vectors upon any axis is equal to the 




algebraic sum of the projections (Fig. 8), whence the proof of the 
distributive law (15). 
Similarly, 

A(B+C+D+E+...)=AB+AC+AD+AE + ... , 
and also 

(A + B)(C+D)=(A+B)C + (A+B)D 

= C(A+B)+D(A+B)=AC + BC+AD+BD. 
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And since B-C is the same thing as B + (-C), we have also 
A(B-C)=AB-AC. 

In fine, the scalar multiplication of vectors is commutative as well 
as distributive, and any two vector polynomials are multiplied out 
precisely as in ordinary algebra. This makes the scalar multiplica- 
tion of vectors a powerful operation. 

As examples we may quote 

(A+B)(A-B)=A a -B a = ^ 2 -5 a , 

meaning that the product of the lengths of the diagonals of a 
parallelogram multiplied by the cosine of their included angle is 




equal to the difference of the squares constructed upon the sides 
of the parallelogram ; again, 

(A+B) 2 = ^ 2 + 5 2 + 2AB, 
or (Fig. 9), remembering that AB =AB cos (x - 6) = -AB cos 0, 

C 2 = A 2 + B 2 -2AB cos 6, 

the well-known trigonometrical relation. In particular, if A J_ B, 

(A+B) 2 = ^ 2 + 5 2 , 

the theorem of Pythagoras. As a third example, let us quote the 
scalar product of two coinitial unit vectors, written as in (7a), 

r x = [i cos fa + i sin fa] sin 8 X + k cos 6 lt 

r 2 = [i cos fa + j sin fa] sin 2 + k c ° s 0g, 

and representing (by their end-points) two places on the Earth * 
whose geographic colatitudes and longitudes are 6 lt fa and 6 2 , fa. 
If s be their geodesic or shortest distance, i.e. the angle (r^ r 2 ), 
we have coss = r x r 2 . Now i 2 =i, etc., and ij=jk = ki=o. Thus, 

* Assumed to be ideally spherical, of radius taken for unit length. 
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multiplying out the two trinomials we have, for the required dis- 
tance s, 

cos 5 = cos 6 1 . cos 2 + sin B x . sin 2 (cos <£ a . cos <£ 2 + sin ^ . sin <i> 2 ). 

Again, calling for the moment a, b two equatorial unit vectors 
having the longitudes of the two places (cf. Fig. 6), viz. 

a = i . cos <f> 1 + j.sin<j> 1 , h=\ cos <j> 2 + j sin <f> 2 , 
we have 

ab = cos (<j> t - <j>^j = cos ^ cos <j> 2 + sin fa sin <j> 2 , 

the well-known formula of plane trigonometry, so that the geodesic 
distance of the two places, (0 X , fa) and (6 2 , fa), becomes 

cos 5 = cos $! cos 9 2 + sin 6 1 sin 6 2 cos (^ x - <j> 2 ), (7b) 

an important formula for navigators, which is at the same time the 
fundamental " cosine formula " of spherical trigonometry. In 
fact, N being the pole (0 = 0), formula (7b) concerns the spherical 
triangle 1N2 (Fig. 10), whose sides are 5, 6 lt 6 2 , and whose angle 

N 




included between the latter two is $ 2 — ^ Notice that this is 
valid for any spherical triangle ; for one of its corners can always 
be considered as our pole, 6 = 0. 

The reader will not be astonished to see the comparatively 
complicated theorems of euclidean geometry thus to follow with- 
out the least trouble from squaring the sum of vectoFsfiir from 
multiplying scalarly two unit vectors. For essentially all euclidean 
relations have been condensed into the above vectorial definitions 
and rules of operations (addition and scalar multiplication). Still, 
as such a condensed system, the vector algebra is exceedingly 
useful. The reader will find for himself that the vector equality 
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and the vector addition alone, as explained in Sections 1 to 4, even 
without the help of the scalar product, are sufficient to demonstrate 
formally a large number of euclidean theorems, such, for instance, 
as the mutual bisection of the diagonals of a parallelogram, the 
common cross of the three medians of a triangle, and so on. 

The scope and purpose of this booklet do not permit us to enter 
into all these attractive details. The willing reader will, however, 
find no difficulty in treating them as exercises which he will soon 
find to be easy as well as interesting and useful, when skill in 
handling the vector method is aimed at. 

6. The vector product of vectors. Two non-collinear vectors, A 
and B, can always be said to define a plane A, B, by making them 
coinitial, for instance, as in Fig. 11. We already know that one 




Fig. 11. 

of the previous operations, AB, deprives them of all their properly 
vectorial ^characteristics, and the other, A+B, or more generally 
xA + yB, gives us only vectors which are again in the plane A, B. 

The operation to be now introduced is in this respect particularly 
interesting, since it yields a vector outside the plane of the operands 
A,B. 

Definition. We call vector product of A into B and denote by 
VAB a third vector C normal to A, B and drawn so that for an 
observer glancing along C the rotation turning A into B, through 
an angle smaller than 180 , is clockwise. This fixes the direction 
and the sense of the vector product C=VAB, and its tensor is 
defined as equal to the area of the parallelogram constructed upon 
A, B as sides, i.e. 

C = |VAB|=AB|sin(A,B)|. (16) 

From this definition we see, first of all, that the vector product 
is not commutative, inasmuch as we have 

VBA=-VAB. (17) 
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Again, if A,B are parallel to one another, i.e. collinear, we have 

VAB=0. 
And if A JL B, then sin (A, B) = ± i, and 
C=\VAB\=AB, 
while A, B and C form a right-handed normal system of three 
vectors. If A points upward and B towards the right, then 
C=VAB points forward. 

If we know of two vectors A, B that their vector product vanishes, 
then we can conclude only that they are parallel (collinear), i.e. 
that 

A = wB, 

where m is some undetermined scalar number, but by no means 
that one of the vectors vanishes (unless we know beforehand that 
they cannot be parallel). This is, mutatis mutandis, analogous to 
what has been said in Section 5 with regard to the scalar product. 
From (16) we see at once that the vector product of m times A 
into n times B is equal to 

wwVAB. 
Thus, for instance, 

VAB = ^BVab, (18) 

where a, b are the units of A, B. Similarly AB =ABaib. 

For a right-handed system of normal unit vectors, as the previous 
i, j, k, we have 

Vij=k, Vjk=i, Vki=j, (a) 

three relations derivable from one another by cyclic permutations 
of i, j, k. At the same time we have, of course, as for every vector, 

Vii=Vjj=Vkk = o. 

Contrast these relations with the previous ones, i 2 =j 2 = k a =i and 
ij=jk = ki = o. The latter follow also from (a) ; for, by the second 
of (a), for instance, i=Vjk is normal to j, and therefore ji=jVjk=0. 
It is scarcely necessary to explain that jVjk means the scalar 
product of the vectors j and Vjk. 

More generally we have, for any two vectors A, B, by the very 
definition of VAB, 

AVAB=BVAB=0. 

Let now A, B, C be any three vectors whatever, generally non 
coplanar with one another. Then the scalarly-vectorial product, 

AVBC, 
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which is itself a scalar, has a very simple geometrical meaning. 
In fact, let A, B, C (in the order as they are written) form a right- 
handed system, i.e. such that a person glancing along C sees the 




rotation from A to B (through less than ir) clockwise. Construct 
upon A, B, C as edges a parallelepipedon (Fig. 12). Then VBC 
will be perpendicular to the base B,C, and its tensor will be equal 
to the area of this base ; in symbols, 

VBC = (area of base) n, 

where n is a unit vector perpendicular to the base. Therefore, 

AVBC = (area of base) An, 

and An being the height of the parallelepipedon, we see that 

AVBC = volume of parallelepipedon A, B, C, 

provided that A, B, C is a right-handed arrangement of the edges. 
(If it were a left-handed arrangement, then AVBC would be equal 
to minus the volume.) Now, the same volume can be expressed 
by taking C, A or A, B as base. Thus we obtain the important 
property 

AVBC=BVCA = CVAB, (19) 

or in words : the cyclic permutation of the three factors of AVBC 
does not influence the value of the product.* Inverting the cyclic 
order is equivalent to changing its sign. For VCB is minus VBC. 
The particular property («A + yBjVAB=0 can now be interpreted 
geometrically by saying that the volume of a parallelepipedon 

* The validity of formula (19) is by no means based upon this volume- 
proof (or rather illustration), which is given here only because it best 
appeals to simple intuitions. In fact, (19) can be proved algebraically, 
without any appeal to the concept of ' volume.' 
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vanishes when its three edges become coplanar, that is to say, 
when all its faces collapse into one plane. 

If of any three vectors A, B, C we know that 

AVBC = o, 
then the only thing we can conclude is that A, B, C are coplanar, 
but by no means that one of these vectors vanishes. Conversely, 
if A, B, C are coplanar, we have AVBC =o. The theorem expressed 
by (19) is of great utility in many applications, and it deserves, 
therefore, to be well kept in mind. 

As in the case of the scalar product, one of the most important 
properties of the vector product is its distributivity, i.e. for any three 

vectors A, B, C, 

VA(B +0) =VAB + VAC. (20) 

This capital property can be proved in a variety of ways. First 
of all, by an immediate geometrical construction of both the 
right- and the left-hand member of (19), — which will be left as 
an exercise for the reader. (It will be enough if the reader con- 
structs it for the simplest case of coplanar A, B, C.*) Another, 
comparatively simple proof, based upon (19), is this: 

Let us write 

VA(B + C)-VAB-VAC=X. 

Then our problem is reduced to proving that X vanishes. Now, 

all the three addends being perpendicular to A, so is their sum 

X, i.e. 

XA=o. 

Again, XB = BVA(B + C) - BVAC 

= (B + C)VBA-CVBA, by (19), 

= BVB A + CVB A - CVB A = o, 

and similarly XC = 0. Thus, the vector X either vanishes or is 
normal to each of the three vectors A, B, C. Now, if these are 
not coplanar, the latter case is excluded, so that X=0. Thus, 
for non-coplanar A, B, C the distributive property (20) is already 
proved. And if A, B, C happen to be coplanar, add to C, for 
instance, a fourth vector D inclined to the plane of A, B, C. Then 
the new vectors A, B, C+D will not be coplanar, and 
VA(C+D)+VB(C+D)=V(A + B)(C + D), 

* For the case of non-coplanar A, B, C is more easily dealt with by the 
following analytical method. 
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and since D can always be so chosen as to make the three relevant 
vectors in each of these products non-coplanar, they may be 
expanded, giving 

VAC+VAD+VBC+VBD=V(A+B)C+V(A + B)D; 

but, by the above, 

VAD+VBD=V(A+B)D, 

whence 

VAC+VBC=V(A+B)C, 

or, changing the sign of both sides, 

VCA+VCB=VC(A+B). 

Thus the distributive property of vector multiplication is proved 
for any A, B, C, coplanar or not. 

The product of two binomials (or polynomials) does not call for 
lengthy explanations. Thus, 

V(A+B)(C+D)=V(A+B)C+V(A+B)D 

= -VC(A+B)-VD(A+B)=VAC+VBC+VAD+VBD. 

The vector multiplication of any two vector polynomials is thus 
seen to obey the same rules as ordinary algebraic multiplication, 
the only difference being that vector products are not commutative. 
A reversal of the order of the two factors changes only the sign 
of their product, which is easily remembered. 

7. Expansion of vector formulae. Basing ourselves upon the 
distributive property just proved, we can at once expand the 
vector product of any two vectors into its cartesian or any other 
form. Thus, if 

A = ^ 1 i + ^ 2 j+^ 3 k, and B=B 1 i + B 2 ]+B 3 k, 

we have, remembering that Vii = 0, Vjk= -Vkj =i, etc., 

VAB=i{A 2 B 3 -A 3 B 2 ) + j{A 3 B 1 -A 1 B 3 ) + k{A 1 B 2 -A 2 B 1 ), (21) 

exhibiting A 2 B 3 -A 3 B 2 , etc. (by cyclic permutation), as the three 
rectangular components of the vector product. Since |VAB| is 
the area of the parallelogram constructed upon A, B as sides, we 
see at the same time that A 2 B 3 -A 3 B 2 , etc., are the areas of the 
projections of this parallelogram upon the planes j, k ; k, i ; i, j, — 
a well-known result which, however, is more easily seen on the 
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vector method. The last formula, (21), is easily memorized in 
its determinantal form, which is 

k 



VAB = 



1 1 

A A A\ 



B* B a 



(2Ifl) 



B\ J->1 *->$ 

In exactly the same way the reader will show himself that the 
cartesian expansion of AVBC, the triple product representing the 
volume of the parallelepipedon A, B, C, is 



AVBC = 
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t-i t-2 t-3 



Bo B* 



(22) 



This, in fact, is the most familiar expression for the volume of 
the parallelepipedon constructed upon A, B, C as edges. Formula 
(22) gives also an immediate verification of the property 

AVBC = BVCA, etc., 
as in (19). For 
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and so on. 

For the scalar product we have immediately, remembering that 
i 2 =;i, ij=o, etc., 

AB = A X B X + A»B„ + A Z B 3 . (23) 

As particular cases of (21) and (23) note the results for two unit 
vectors a, b which include the angle sr, 

sin 2 © = (0^3 - fla&g) 2 + (a^ - o^) 2 + (ajb a - o^) 2 . 

COS TS = Ojbi + fl 2 6 2 + fl S^3> 

«!, b lt etc., being now the direction-cosines of a, b relatively to 
i, j, k as axes. For such is the meaning of the components of 
unit vectors. 

In order to give at least one illustration of the utility of AVBC, 
let us consider three coinitial unit vectors whose end-points may 
be conceived as the vertices 1, 2, 3 of a spherical triangle drawn 
on a unit sphere. Let us use the colatitude and the longitude as 
in (7a). Without any loss to generality we may put the pole 
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(0 = 0) into the vertex I and take the first meridian along the 
side 12 ; thus, 04 being the angle at I, and s 2 , s 3 the sides of the 
spherical triangle opposite 2 and 3, 

ri = i, 

r 2 =i cos s 3 +j sin ,? 3 , 

r 3 =icos 5 2 + sin s 2 . [j cos a x +k sin a-J. 

This gives for the scalarly-vectorial product, by (22), since the 
first vector has no second and no third component, and the second 
vector no third component, 

TjVtgz = sin s 2 sin s 3 sin a Xl 

which, by the cyclical property (19), is also equal to r 2 Vr 3 r a and 
to r 3 Vr 1 r 1! , and these products are obviously equal to 

sin s 3 sin s t sin o 2 
and to 

sin s t sin s 2 sin a 3 , 

where oj, a 3 are the remaining two angles of the spherical triangle. 
Thus, 

sin (4 ^ sin a^sing,, ,j- » 

sin 5 X sin s 2 sin s 3 ' 

the fundamental " sine formula " of spherical trigonometry, 
following on the vector nrethod as easily as the " cosine formula " 
given before. It is interesting to note that the " sine formula " 
is, in this circle of ideas, but the statement of the triple expressi- 
bility of the volume of the parallelepipedon t v r 2 , r 3 , viz. as iJJtga, 
or TjjVrgT! or rgVrjrj. Other examples are left to the care of the 
reader. 

8. Iteration of vectorial multiplication. There is but one more 
important formula to be noted in connection with the vector 
product of vectors, viz. a formula giving a convenient vector 
expansion of the result of repeated or iterated vector multiplication, 

VA(VBC) or simply VAVBC, 

which reads : having obtained the vector, product of B, C, multiply 
it, again vectorially, by A. This ternary product, which occurs 
very often, is, of course, again a vector, to wit, perpendicular to A 
and to VBC ; but the latter being itself perpendicular to B, C, 
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our new vector VAVBC is coplanar with B, C, so that we know 
beforehand that the result will be of the form * 

VAVBC = £B+yC, 

where /8, y are some scalars. Since the ternary product is per- 
pendicular to A, we have /8(AB) +y(AC) =0, so that 

VAVBC = A.{B(CA) -C(AB)}, 

where A. is a scalar. It remains to determine its numerical value. 
This can be done, for instance, in the following manner. First 
of all, A can always be assumed to be coplanar with B, C, since 
its part normal to B, C contributes nothing. Next, dividing both 
sides by ABC, the equation becomes 

VaVbc = A.{b(ca)-c(ab)}, 

where a, b, C are the units of A, B, C. Now, multiply both sides 
scalarly by b, and notice that, by (19), 

bVaVbc = (Vbc) (Vba) =sin (b, c) . sin (b, a). 
Thus, 

sin (b, c) . sin (b, a) = A.;[cos (c, a) - cos (a, b) . cos (b, c)] ; 

but, the three vectors being coplanar, we have 

cos (c, a) =cos (b, a) . cos (b, c) +sin (b, a) . sin (b, c), 

so that A. = i. 

The required formula is, therefore, 

VAVBC = B(CA)-C(AB). (24) 

As an exercise, the reader may verify it by an iterated application 
of the cartesian expansion (21) or (21a). Having once obtained 
this important formula, there will be no difficulty in dealing with 
quaternary vector products, as VDVAVBC, which becomes 
(CA)VDB-(AB)VDC, etc. But such products will hardly occur 
in practice. 

A notable property of the above ternary product and of its two 
cyclical permutations is that 

VAVBC + VBVCA + VCVAB = o, (24a) 

identically. For the six right-hand terms of (24) and of the two 
similar equations destroy themselves in pairs. 

* The trivial case of B, C collinear can be discarded ; for then VAVBC =0. 
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A particular case of (24] which often occurs is that in which C 
is equal to A and is an unit vector u, say. Then we have 

VuVBu = B-(Bu)u, {24b) 

whence we see also that VuVBu is the part of the vector B normal 
to u, in both size and direction. For (Bu) u is the part of B along u. 

To close this section, and at the same time the essential part 
of the whole Vector Algebra, but a few more remarks which will 
be useful in connection with problems often occurring in practice. 

Let X be an unknown, and A, u two given vectors, the latter 
an unit vector. If we know of X only that 

VXu = A, {a) 

we cannot fully determine X. For to a solution of this equation 
we can add any vector rau (since Vuu = o), and X + mu will again 
be a solution of this equation. In order to determine X uniquely 
we must have one more (scalar) datum. Let this be 

Xu=o-, (b) 

where o- is a given scalar. Then X is completely determined. In 
order to find its value explicitly in terms of the given A, u, <r, 
multiply the equation (a) vectorially by u ; then, in virtue of (24), 

X-(Xu)u=VuA, 
and by (b), 

X = <ro+VuA, (c) 

which is the required solution. This simple rule, (c), for solving 
the equations (a) and (b), will often be found helpful. 

9. The Linear Vector Operator. Let R be a variable vector, 
that is to say, one that can assume in turn all possible sizes (tensors) 
and directions. Of each of these determined vectors we can speak 
as of the special value of the variable R. To have a good picture 
of such an abstract concept, imagine R as a straight, extensible 
and contractile string fixed at one of its ends at a permanent 
point ; then its free end-point P occupying in succession all 
possible points of space, OP will represent the various values of R. 
The vector R can, in such a connection, be advantageously called 
the position vector of the point or, if one prefers, of the particle P. 
Now imagine that there is another particle P', and let its position 
vector, with the same origin 0, be called R'. Let there be some 
mechanism, or else our own imagination, which to every chosen 
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position of P makes correspond a certain position of P'. This 
we may express by saying that to every value of R corresponds a 
certain value of R', by writing 

R'=5TR, (25) 

and by calling R' a vector function of the variable vector R. If, 
as we assume, to every R corresponds but one R', determined in 
size and direction, we will say that R' is a monovalent function 
of R, and we will call XS a monovalent vector operator, the symbol 
of some operations to be performed on R in order to obtain R'. 
We can think of such operations in the algebraical, as well as in 
the physical sense of the word, as turning round the representative 
string, stretching or contracting it according to some more or less 
complicated prescription. It is needless to explain that an equation 
such as (25) is equivalent to three scalar equations : each of the 
components of R' equal to some function of, in general, all the 
three components of R. 

Suppose now that R is represented as the sum A+B of some 
two vectors. In general the operations embodied in XS may be 
such that ct(A+B) is not the same thing as ctA+CTB. A good 
example of such an operator is that which converts an incident 
luminous ray into the refracted ray (cf. Simplified Method, quoted 
in Preface). But the operations represented by XS may also, in 
particular, be such that 

xs{A.+B)=xsA. + X3B, 

whatever the vectors A and B. If such be the case we call XS a 
distributive operator. An example of this kind is afforded by£he 
" reflector," i.e. that operator which converts the incident ray into 
the reflected one. The simplest example of a distributive vector 
operator is, however, a scalar number <r used .as a factor ; for we 
have, of course, 

cr(A+B)=crA+<rB. 

This operator is a pure stretcher or (if |<r| < 1) a contractor, and, 
if <r < o, an invertor at the same time. 

Leaving these examples, let us turn to the general distributive 
operator, of which we will only assume that it is a continuous 
operator, i.e. that XSR is a continuous vector function of R. Such 
distributive operators have very far-reaching applications in 
many branches of geometry and physics. They are known better 
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under the name of linear vector operators, and were first introduced 
by the great Hamilton. In fact, it can easily be shown that the 
continuous and distributive operator O when applied to a vector 
R yields another vector R', whose components are each a linear 
function of the components of R, whatever the triad of axes em- 
ployed for the decomposition. For, if n be any integer positive 
scalar number, we have, in virtue of the assumed distributive 
property, 

SJ(mA)=wct(A), 

and this property can easily be extended to negative and fractional 
values of n, and ultimately, by the often repeated limit-reasoning, 
to any real value of n. This being granted, let a, b, C be any 
non-eoplanar unit vectors ; then, whatever R, we can represent it by 

R = aa+yb+zc, 

where x, y, z are Some scalars, the components of R. Thus, the 
corresponding R' will be 

R' = 7S(xtL + yb + ze)=xT3&+ytsb+ziSe, (26) 

and, therefore, the components of R', i.e. #'=aR', etc., 

#' = (aETa)# + (aE7b)y + (asJc)z, ) , , k 

y (26a) 

y' = (bwa)* + etc. ; 3' = (ecra)# + etc.,J 

and these are linear, homogeneous functions of x, y, z, the com- 
ponents of R, — the coefficients (aCJa), etc., being certain scalars 
(scalar products of a and GTa, a and wb, and so on), which depend 
partly on the nature of the operator TS, and partly on the choice 
of the framework a, b, c. 

From (26a) we see that the operator TS is fully determined if 
we give, for any chosen a, b, C, the nine coefficients 

ana, acrb, anc, boa, etc., 

which we will denote by 

&aa, EU> SToc. w JO , etc., 

respectively. And it is not difficult to see that these nine scalar 
data, which are sufficient, are also necessary to determine com- 
pletely a linear vector operator. 

We can express the same thing more simply by taking the vector 
equation (26) instead of its components (26a). Let us rewrite (26), 
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taking for a, b, c any triad of normal vectors ; then x, y, z stand 
for aR, etc., so that the equation is 

R'=wa(aR)+crb(bR)+crc(cR). (26O 

Thus we can say that the operator © is fully determined if we 
know what it yields if applied to the three conventional vectors 
a, b, c, that is to say, if we are given the three vectors 

era, crb, ere. 

Each of these implies 3 scalar data, so that in all we have again 
9 data, as before. What we have, a moment ago, denoted by 
ts^, CT oS , etc., are simply the cartesian components of these three 
vector data, taken along the axes a, b, C. 

Hitherto we have spoken of the most general linear vector 
operator. Let us now explain an important subdivision of this 
vast class of operators. Let A and B be any two vectors. We 
can take ctB and multiply it scalarly by A, or first form JtrA and 
then multiply it so by B. In this way we should obtain the two 
scalars, 

AwB and BstA. 

Now, in general, these two numbers will be different from one 
another. But the operator xs may happen to be such that they 
are equal, i.e. that, for any A, B, 

AWB=B©A. (27) 

If such be the case we call ST a self-conjugate operator. By this 
definition (27) we shall also have VS ca = ZSi a , etc., so that a self- 
conjugate operator has but six mutually independent (that is, six 
independently prescribable) coefficients, or constituent's, 

Goo. ™ab ™ „c 

rs K . 

This table, after the insertion of xs ba = rs ab , etc., at the vacant 
places, is symmetrical with respect to its diagonal ; whence also 
the name of symmetrical operator, used as a synonym for the self- 
conjugate operator. 

This being a sub-class of TS, the general operator, the remainder 
of the class of C7's, for which 

E^Wj.. etc., 
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are called non-symmetrical or asymmetrical operators. If vs be 
such an operator, and if vs' be another operator such that 

AsjB=Bst'A, (28) 

for any pair A, B of vectors, then vs' is called the conjugate of VS. 
Obviously also VS is the conjugate of VS'. Thus, if vf^, etc., be 
the coefficients of VS' ', we have, remembering that VS^ stands for 
acrb, etc., 

^06 = ©fa,, ^'6 =W o6 , CT' co =W ac , (28a) 

while, of course, V5 , aa =T3 <m , etc. Thus we see also that to every 
operator VS there is one (and only one) conjugate vs'. In particular, 
if © is a symmetrical operator, its conjugate is identical with it, 
whence " self-conjugate " as a synonym of symmetrical operators. 
In harmony with this, (27) is but a special case of (28), viz. for 
VS' = VS. 

Let us use vs for any linear vector operator, and to for symmetrical 
operators only. (In fact, without the circumflex this last letter 
of the Greek alphabet has some symmetry.) 

Manifestly the symmetrical operator &> will be a great deal 
simpler than the asymmetrical VS. It is, therefore, very agreeable 
to see that any VS can be split into an o> and some other asymmetrical, 
but very simple, operator which is called an antisymmetrical (or 
skew-) operator and which we will denote by a. The latter is 
denned most conveniently by saying that, for any A, B, 

AaB=-BoA, .'. AaA=0, (29) 

and therefore also a ai = -o te , etc., and a M = o, etc., so that the 
table for such an operator becomes 

-"■at a bc (29a) 

which justifies the name. The announced property can shortly be 
written 

VS = a) + o, 

which is a symbolic short for 

CTR = a)R + aR, 

where R is any vector operand. The said property is easily 
proved. 
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In fact, let 75' be the conjugate of the given operator ST. Then 
we have, identically, 

VS^TS + TS^+^TS-TS'). (30} 

But the first term represents a symmetrical operator, because, 
by (28), 

A(CT + cr')B = AwB + Asi'B =Bcr'A +BCTA =B(ST + ST')A, 
which precisely is the definition (27) of a symmetric operator. 
And the second term is antisymmetric, for 

A(st-gt')B=Bct'A-BstA = -B{a-rs')A, 

as in (29), the definition of antisymmetric operators. This proves 
the statement, without the slightest need of splitting ZS into its 
nine constituents rs m , etc. 

We thus see that every linear vector operator can be written 

ST = co + o, (31) 

where its symmetrical part is eo = J(gt + gx') and its antisymmetrical 
part o = |(ct-ct'). 

If the reader so desires he can introduce the nine coefficients of 
these operators. Then 

proving again that <o is self -conjugate, and 

a aa = 0, etc. , a ai = 1(13^ - ©J = - o^, 
proving that a is antisymmetric. 

Turning now to the antisymmetric operator a we can see from 
its definition (29) that it has a very simple meaning. Jn fact, 
let R be the vector operated upon. Then, by the second of (29), 
whatever the value of R, aR is a vector normal to R. Now, this 
condition can be satisfied by putting 

aR=VwR, 
where w is some fixed vector. But such being the case, we have 
also, for any A, B, 

AaB = A VwB = - BVwA = - BaA, 

so that the general definition (29) is completely satisfied. Thus, 
the antisymmetric operator is, dropping the arbitrary operand, 

o=Vw ; 
in words, to operate with a is to multiply vectorially by a certain 
vector w. 



THE LINEAR VECTOR OPERATOR 31 

Ultimately, therefore, we can write, instead of (31), for any- 
linear vector operator, 

CT = w + Vw; (32) 

the symmetric operator <o is half the sum of TS and of its conjugate 
ST', while w is a certain vector characterizing the given operator TS. 
Notice that w, being symmetric, implies 6 independent scalar data, 
and w, being an ordinary vector, 3 more, making in all 9, as before. 
Obviously, o> and w can be prescribed independently of one another, 
and these two data (equivalent to 6 + 3=9 scalar ones) fully 
determine the asymmetric operator TS. 

If we desire to express w in terms of the coefficients TS ai , etc., 
we can easily do so. For, from the table or the " matrix " (29a) 
we see that ' 

oR =&{a ai R b + aji c ) + etc., 

and since <* (l4 = |(sT lli -B7j a ), and so on, while aR = VwR, we find 
without difficulty that, if a, b, c be a ng^-handed system, 

2w=a(sr c6 -w 5c ) +\>{TS ac -TS ca ) +o(sx i .-wJ, (33) 

which is the required expansion of w. 

Having thus shown that the antisymmetric part of any operator 
TS is simply a vectorial multiplier Vw, it will henceforth be enough 
to study the remaining part of TS, that is to say, the symmetrical 
operator co. 

Principal axes of 01. — Let R be the operand. Then the vector 
R' =(oR will in general differ from R not only in size but in direction 
as well. But if R assumes certain particular directions, then it 
may happen that R' coincides with R in direction, if not in size. 
Let x, an unit vector, represent such a privileged direction. Then, 
(oR being a linear function, the inverse direction - x will, obviously, 
partake of the same privilege. Such particular directions ±x are 
called principal axes of the symmetrical operator co, both +x and 
-x counting for one axis. 

This is merely a definition. Let us now see whether at all such 
axes and how many of them do exist, and what are their mutual 
relations. Let us start with the last question. Suppose then that 
there are two different principal axes x and y. Then, by the very 
definition of such axes, 

WX=0,1 H (34) 
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where Wj, a> 2 are some ordinary scalar numbers, which are called 
the principal values of u, corresponding to these axes x, y. Multiply 
the first equation scalarly by y, and the second by x, and subtract 
them from one another. The result will be 

(<!>■! - <o 2 )xy = y">x - x<uy. 

But, the operator m being symmetrical, xa)y=y<i>x. Thus 

(a) 1 -(o 2 )xy = 0, 

and if o^ + co 2 , we have xy = 0, that is to say, x _L y. And should 
it happen that ojj^coj, i.e. 

(ox = (OjX, <uy = a^y, 

then X, fi being any two scalar numbers, 

t>)(A.x + /iy) = w^Ax + py). 

But Ax+/*y is any vector in the plane x, y. We thus see that if 
there are two principal axes x, y to which correspond different 
principal values o^, «u 2 , these axes must be normal to one another, 
And if <o 1 = ui 2 , then every direction in the plane x, y is also a 
principal axis. 

Suppose now there is still a third principal axis z not coplanar 
with x, y, and let g> 3 be its corresponding principal value, so that 

b)Z = OJfjZ. 

Then, reasoning as before, we shall see that if o) 3 , u> lt <a 2 are all 
different, z will be normal to x and also to y. And if gij = oj 2 = a> 3 , 
then every direction whatever will be a principal axis with the 
same principal value, in which case the operator <o degenerates 
into an ordinary scalar factor. 

Thus, in the most general case the symmetrical operator to can 
have three different,* mutually perpendicular principal axes, x, y, 
z ; and only three. Because the fourth, if it existed and carried 
a new a> 4 , would have to be normal to those three which, in our 
space, is nonsense ; and if a> 4 were equal to o^, say, then the whole 
plane passing through the fourth and the first axis would consist 
of principal axes, and since this plane would cut the y, z plane, 
ci> 2 and (i) 3 could not be different from one another, against the 
assumption. 

* I.e. such to which correspond different principal values. 
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Having thus settled the question about the number of the 
possible different principal axes of a> and their mutual orientation, 
it remains to see whether they exist, or better, to find them. The 
technical side of the latter problem will depend upon the manner 
how a) is given. Suppose it is given through its six different 
coefficients 

<"<», <»bb, <"«; "<*, »>tc, Ka 

with respect to some arbitrarily fixed framework of normal unit 
vectors a, b, C, or — which is the same thing — that the three vectors 

o>a, wb, we 

are given, say, equal to A, B, C, respectively, so that (co being 
symmetrical) Ab = Ba, etc. Let x be a principal axis and n the 
corresponding principal value (both to be found). Then if x x , x 2 , 
x 3 are the direction cosines of x with respect to a, b, c, so that 

x=x 1 a,+xjb+x 3 c, 
we have 

ox = x x u>% + x 2 o)b + x 3 o)C = x x h + x^i + x 3 C, 

and since a>x = wx, 

%A + x^i +x 3 C = n (x-fi + xjb + x 3 c) 
or 

%(A -wa) +x 2 (B-nh) +x 3 [C-nc) =0. (35) 

From this equation we see that the three vectors A-na, etc., are 
coplanar, so that the volume of the parallelepipedon constructed 
upon them is nil, i.e. 

(A-rca)V(B-«b)(C-«c)=o. (36) 

Since A, B, C are given, this is a cubic equation for the unknown n. 
Multiply it out and remember that a = Vbc, and therefore aVbc = 1 . 
Then the result will be 

n 3 - w 2 (Aa + Bb + Cc) + w(aVBC + bVCA + cVAB) - AVBC = 0. (36a) 
Each of the coefficients of this cubic equation for the principal 
values of the operator to has a simple geometric meaning : the 
first is the sum of the projections of the vectors A = wa, etc., upon 
the conventional a, b, c, the second the sum of the volumes of 
the parallepipeda a, B, C, etc., and the last is the volume of the 
parallelepipedon A, B, C. At the same time we see that these 
three expressions are invariants of o>, i.e. independent of the choice 
of the reference system a, b, c. In fact, if %, w 2 , n 3 be the principal 

S.V.A. c 
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values of a>, which manifestly are intrinsic properties of the operator, 
independent of the reference framework, we have, by (36a), 

Aa+Bb + Cc = w 1 + w 2 +w 3 , 
aVBC + bVCA+cVAB=n 2 W3 + M 3 M 1 + m 1 »j 2 , - (37) 

AVBC = n 1 n i n 3 , where A = <ua, etc. 

These are very important formulae, exhibiting the three invariants 
of the symmetrical operator o. 

Now, if only A, B, C are real, as we assume, all these invariants, 
i.e. the coefficients of the cubic (36a) are real. That equation has, 
therefore, at least one real root. Let this be Mj, and let us take 
the corresponding principal axis * as our reference axis a. Then 

A = ua=M 1 a; .'. bVCA=w 1 Ce, cVAB = w 1 Bb, 
and the left-hand member of (36a) becomes at once 

n 3 - nh^ + (ntt! - n 2 ) (Bb + Cc) + (n - ra^aVBC, 

which is, as it should be, divisible by n - n^, leaving for the remain- 
ing two principal values w 2 , n 3 the quadratic 

n* - «(Bb + Cc) + aVBC = 0, 
which gives 

"» = I (Bb + Cc) ± s/i(Bb+Cc) a -aVBC, (38a) 

71 3 

or, in terms of the coefficients u u = bB, etc., since 



£*=l(«to+«0 ±>/iK-<0 2 +<, (38) 



so that, if only all the coefficients <■>„ are real, these two principal 
values and, therefore, also the corresponding principal axes are 
real. That they form with the first axis a normal system we 
already know. 

We have written down the two roots (38) in the assumption 
that a) was given by prescribing its coefficients v> M or the vectors 
A, B, C, with respect to an arbitrary framework a, b, c. But, 
as a matter of fact, this expansion of the roots is superfluous. 
For, having taken a as one of the principal axes of a>, we know 
beforehand that b, C will be its remaining two axes, i.e. that 
B = o)b = M 2 b, and C=w 3 c. 

* Whose direction cosines with respect to any a, b, c might at once be 
determined from (35) by taking in it n =n v 
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Now, with these values, we have aVBC = w 2 w 3 aVbc = w 2 w 3 , so that 
(38a) becomes 

5 = l(«i + »»)±Vi(n,-« I )« 

which is, as it should be, an identity. Thus, the only necessary 
thing was to state that the cubic (36a) has at least one real root, 
and this was immediately clear. 

Having thus ascertained the general properties of the principal 
axes of to, let us take them as our (natural) reference system 
a, b, c, which we will now call i, j, k. Then, n±, n^, n s being the 
corresponding principal values, the most general symmetrical linear 
vector function will be 

«R = w 1 i(iR) +WjjJ(J R ) + w 3 k(kR), 
that is, n± times the first component of R along i plus, etc., or 
using the dot, instead of brackets, as separator, 

<oR — n^. . iR +n 2 j . jR + « 3 k . kR, 

or dropping the operand R, which it is useless to repeat so many 
times, 

w = w 1 i.i + M 2 j.j+w 3 k.k. (39) 

Thus the symmetrical operator assumes the form of what is called 
(after Gibbs) a dyadic, which is a polynomial, in our case a trinomial, 
of dyads such as n^.i, etc. It will be well to say a few words 
on these useful mathematical beings. 

Dyads and Dyadics. — The dyads appearing in (39), which, apart 
of the scalar factors (calling for no explanations), are of the form 
i . i, are but special cases of the general dyad which is 

a. b, 

a, b being any two vectors, in general that is not coinciding in 
direction ; the first vector is called the antecedent, and the second, 
the consequent of the dyad. The dyad as an operator can be used 
either as a prefactor of the operand, say 

a . bR, meaning a(bR), 
or also as a postf actor, 

Ra . b, which means (Ra) b. 
As Heaviside says somewhere : "A cart may either be pulled or 
pushed." In fact, this two-fold possibility of attelage of the 
operator turns out to be very advantageous. 
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If a, b are not collinear, then a . bR is, of course, altogether 
different from Ba . b. Such is the case with the most general (asym- 
metric) dyad. But if the antecedent and consequent happen to 
be collinear, as in the case of (39), then the dyad, applied to any 
vector, yields the same result whether it acts as a pre- or a post- 
factor. Such dyads are called symmetrical dyads. They are all 
of the form 

era . a, 

where a is an unit vector, and o- a scalar. Sums of such dyads 
are called symmetrical dyadics. Thus, the most general sym- 
metrical or self-conjugate linear vector operator <o may be repre- 
sented as a (trinomial) symmetrical dyadic. Such, in fact, is (39). 
Consider any, generally asymmetric, dyadic, say a trinomial one, 
<£=a.x + b.y+c.z, 

which is a certain linear operator. Interchange the antecedents 
and the consequents ; then the resulting operator or dyadic 

<£'=x.a+y .b+z.c 
is called the conjugate of <f>. For such it is according to the previous 
definition. In fact, if R, S be any two vectors, we have 

Ra.xS=Sx.aR, etc., 

since the scalar product is commutative. Thus also the product 
of R into <f>& is seen to be identical with that of S into <£'R, — as 
in the definition (28). Again the product of R<£ into S is identical 
with that of R into <£S. Thus no brackets or other separators 
are needed, and the last property can be written simply 

R$5 = S<£'R, 

valid for any dyadic <£, and its conjugate <j>'. 

We have already seen that the self-conjugate linear operator 
<o can be represented as a symmetrical dyadic. We may still 
mention that the general linear vector operator sr can always be 
reduced to what is called a normal (trinomial) dyadic, i.e. 

ST=o- 1 l.i+o- 2 m. j -fo-jn.k, (40) 

where cr 1 , cr 2 , <r 3 are scalars (either all positive or all negative), and 
both the antecedents 1, m, n and the consequents i, j, k form 
normal, say right-handed, systems of unit vectors. If these are 
distinct from one another, we have an asymmetric operator XS 
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and if they coincide, we have a symmetric operator m, as before. 
The conjugate of the general T3 will be 

W = «r r i . 1 + ctj4 .m+o-gk.n. (41) 

The special symmetric dyadic 

i = i.i + j. j+k.k 

leaves, of course, any vector operand R intact, and is, therefore, 
called an idemf actor. It is also, for all purposes, equivalent to I. 
And if a- be any scalar, then o-i as an operator is equivalent to cr 
itself, as an ordinary numerical factor. Thus, expressions such 
as <r + a . b will again be dyadics, and require no further explanations. 

To close this section it will be enough to make a few remarks 
on the "multiplication," i.e. the successive application of dyadics. 

If $ = a . b and ^=c . d be two dyads, and R any vector operand, 
we have obviously 

*(*B) = (W)R, 
where <f>\f/ is the dyad a (be) .d = (bc)a.d. Similarly, if y be a 
third dyad, we have 

W(yB)=*(ftOB-(W)yR, 

the associative property, so that each of these expressions can be 
simply written <f>\f/yR. And the same is easily seen to hold if 
<f>, y\>, etc., stand for binomial or polynomial dyadics. Again, 
since the scalar product of vectors is distributive, we have for 
any dyadic <f> and any vectors R, S, 

<£(R+S)=<£R + #3, 
and also, if xf/, y be two more dyadics, the operational equation 

4>{^ + y)=<f>tf' + 4>y, 
and also 

(^ + -y)<£ = #+y<£. 

In short, the distributive property holds for the multiplication of 
any polynomials of dyads, and therefore of dyadics. Such products 
can, therefore, be expanded as in ordinary algebra, the only neces- 
sary precaution being to keep the order of the operators and of 
the constituents of the dyads intact, since (in general) the com- 
mutative property does not hold. Thus, for instance, 

(a . b + c . d) (e . f +g . b.) = (be) a . f + (bg)a . h + (de)c . f + (dg) e . h. 
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Vectors not separated by dots are fused into scalar products, as 
(be), (bg), etc., and here of course the order is irrelevant ; but it 
must be carefully preserved in the resulting dyads, such as a . f, 
not f.a (unless a, f are collinear). Apart from this precaution, 
the multiplication of dyadics is as easy and convenient as the 
common multiplication of polynomials, and it will be found to 
render inestimable services in the treatment of many geometrical 
and physical, especially optical, problems. Some illustrations of 
the latter kind will be found in the " Simplified Method, etc.," 
mentioned before. The final result of such multiplications of two 
or more polynomials will be a polynomial of dyads, say 

A.B+C.D+E.F+G.H + etc. ; 

but since each of these antecedents and consequents can be 
expressed in the form x& + yb +zc, where a, b, C are any non-coplanar 
unit vectors, any such result can, in the first place, be reduced to 
a sum of nine dyads, viz. 

<r u a . a + o-^b . b + o- 33 c . c + o- 23 b . c + <r 32 c . b + . . . + c^b . a, 

and it can be proved that this can always be reduced by a proper 
choice of two orthogonal systems, i, j, k and I, m, n, to the normal 
form (40), which is that of any, generally asymmetric, linear 
operator rs. Ultimately, the latter can with advantage be split 
into a symmetric operator o> and the simple operator Vw, as 
in (32). 

10. Hints on Differentiation of Vectors. The concepts of 
differentiation and integration as applied to vectors do not 
belong to the subject proper of this booklet, which is vector 
Algebra. Yet a few elementary remarks on the differentiation 
of vectorial expressions may be here added, as they are likely to 
be useful to some readers, and as they do by no means require 
much space. 

Let R be a variable vector. To have a possibly desirable picture 
think of R as the position-vector of a particle moving about in 
space, round a fixed origin. Let t be any independent scalar 
variable, say the time. Then, AR being the vector increment, 
i.e. the vector drawn from the position of the particle at the instant 
t to that at a later date t + Ai, the quotient AR/Ai will be a certain 
vector, having a definite tensor (size) and a definite direction. 
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We may call it provisionally the average vector-velocity of 
the particle. If this quotient (a vector) tends, with indefi- 
nitely decreasing At, to some definite limit, definite both in size 
and direction, we call this limit-vector the derivative or the 
fluxion of R with respect to t (or the vector velocity of the 

particle), and denote it by -j- or R. In short symbols 

dR ■ T . AR 

lF =R=Llm AT 

This vector will, in our illustration, be tangential to the orbit of 
the particle, and its tensor will represent the particle's speed. 
From this definition it follows at once that 

where R, S are any vector functions of the variable t. And, if 
t be the unit of R, so that R = i?r, we have of course 

Again, since the scalar product of two vectors is distributive, so 
that A(RS) =RAS+SAR plus terms of higher -order, we have 

j(RS)=RS+RS. 

In particular, if r be a unit vector, so that r 2 =i, we have, by 
differentiating the latter condition, rr = 0, so that r _L r, which is 
also an obvious property. Similarly, for the vector product, which 
again is distributive, 

^VRS = VRS + VRS, 

the only precaution being to preserve the order of the factors, or 
— if this be inverted — to change the sign of the product in question. 
In quite the same way we have 

~ AVBC = AVBC + AVBC + AVBC, 

at 

and so forth. 

Even the case of linear vector functions, such as C7R, does not 
call for lengthy explanations. If not only the operand R, but 
also the nature of the operator ZS varies with t, we have 

-5-(nrR) = £rR + crR, 
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since ST is distributive. Here rs is the derivative of the operator. 
If, for instance, zs is represented as a dyadic, say 

A.B + C.D + E.F, 

we have 

ct = A. B+A.B + ...+E.F. 

And if the form o>+Vw is used, we have 

W = o> + Vw, 

where o> can again be expanded, as the derivative of a symmetrical 
dyadic. It is scarcely necessary to add any further explanations. 
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